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Spherically symmetric, null dust clouds, like their time-like counterparts, may collapse classically 
into black holes or naked singularities depending on their initial conditions. We consider the Hamil¬ 
tonian dynamics of the collapse of an arbitrary distribution of null dust, expressed in terms of the 
physical radius, R, the null coordinates, V for a collapsing cloud or U for an expanding cloud, the 
mass function, m, of the null matter, and their conjugate momenta. This description is obtained 
from the ADM description by a Kuchar-type canonical transformation. The constraints are linear 
in the canonical momenta and Dirac’s constraint quantization program is implemented. Explicit 
solutions to the constraints are obtained for both expanding and contracting null dust clouds with 
arbitrary mass functions. 
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I. INTRODUCTION 


Spherically symmetric dust clouds, depending on their 
initial matter and velocity distributions, will collapse in 
classical general relativity to form either black holes or 
naked singularities. Black holes are better understood 
than naked singularities. They are generally expected to 
evaporate via their associated Hawking radiation [jl|], al¬ 
though no agreement has yet been achieved regarding the 
end state of collapse, i.e., whether a remnant survives or 
whether all the matter contained in the original cloud is 
thermally radiated away. If a portion of the collapsing 
matter does manage to form a stable black hole, it is ex¬ 
pected that the total mass of the remnant will be quan¬ 
tized. On the other hand, if all the matter is radiated 
away before a stable end state can form then one must 
explain what happens to the information that was con¬ 
tained in the initial matter distribution. The formation 
of black holes therefore presents a number of deep puz¬ 
zles and various approaches to quantum gravity are being 
employed to address these at the present time [|| ||, [| |) . 
On the contrary, naked singularities have received com¬ 
paratively little attention. Yet, the formation of naked 
singularities (singularities that are visible either locally or 
asymptotically) is much more difficult to understand and 
for an entirely different reason: their existence implies the 
absence of a well defined Cauchy problem to the future 
of some light-like surface (the Cauchy horizon), therefore 
any attempt to describe the system to the future of this 
surface fails for lack of initial conditions. It seems that 
space-time must be terminated at the Cauchy horizon. 
In order to avoid the associated problems, Penrose pro¬ 
posed a Cosmic Censor ||, whose function is essentially 
to ensure that naked singularities never form. The mech¬ 


anism by which the Cosmic Censor operates, however, is 
still shrouded in mystery. The Censor is most likely not 
classical because most models of classical collapse lead to 
the formation of both black holes and naked singularities 
in different domains of the initial phase space [7j] . In fact 
very little is currently understood about the final stages 
of a collapse that leads to the formation of a classical 
naked singularity. 

There are indications from the semi-classical treatment 
of naked singularities, in which the gravitational degrees 
of freedom are considered to be classical, that Penrose’s 
Cosmic Censor may, in fact, be the quantum theory it¬ 
self ||. However, at the very final stages of collapse it is 
not possible to treat the gravitational degrees of freedom 
classically and a full blown quantum theory of the grav¬ 
itational field becomes necessary to establish this possi¬ 
bility U firmly. Singularities in general relativity signal 
a breakdown of the classical theory, a regime in which 
the classical equations are meaningless. Cosmic Censor¬ 
ship probably points to the need for quantum gravity in 
the same way as, more than eighty years ago, the elec¬ 
trodynamic instability of atoms pointed to the need for 
quantum mechanics. A good question is just how com¬ 
plete a theory of quantum gravity is required to begin 
addressing such issues as the Cosmic Censor. We take 
the attitude that, from past experience, it is not unrea¬ 
sonable to expect many of the key effects of quantum 
gravity to be understood from a more naive quantiza¬ 
tion of the gravitational field which, while it may be in¬ 
complete, incorporates the essential features of quantum 
mechanics. 

This is what we propose to do in this paper. Our ob¬ 
jective is to consider the midi-superspace quantization 
of a spherically symmetric cloud of null matter speci- 
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fied by an arbitrary mass distribution and collapsing in 
its own gravitational field. The model we are concerned 
with therefore is a solution of Einstein’s equations with 
pressureless, null dust 0 described by the stress energy 
tensor T = e(x)U^U v , where e(x) is the energy density 
of the cloud and U 2 = 0. When the cloud is contracting, 
the solution is characterized by an arbitrary function, 
miy ), of the advanced null coordinate, V £ (— oo,+oo). 
The mass function is generally taken to be vanishing for 
V <V a and constant, M, when V > Vi. The space-time 
is described by the metric 

ds 2 = ^1 - 2TV 'j! > 'j dV 2 - 2 dRdV - R 2 dn 2 , (1.1) 

where R £ [0, +oo) is the area radius. The region V >V\ 
is a part of the Schwarzschild space-time. In this region 
the metric may be written in terms of the Eddington- 
Finkelstein coordinates, U and V as 

1 _ if) 2 ’ ( L2 ) 

where U = V — 2 R* and R* is the tortoise coordinate. 
The region V < V 0 is a part of Minkowski space-time, 
with metric 



ds 2 = dudv — R 2 dfl 2 , (1.3) 

where u and v are the ordinary retarded and advanced 
times, respectively T R. 

In the time reversed situation the null dust cloud is ex¬ 
panding instead of contracting and the solution is written 
in an analogous fashion, 

ds 2 = (\ - dU 2 + 2dRdU - R 2 dn 2 (1.4) 

in terms of a retarded null coordinate, U. Again, the 
mass function is generally taken to vanish when U > 
Ui, having some constant value, M before some earlier 
retarded time, U < U 0 . The region U > U\ is then a 
part of Minkowski space-time, while the region U < U 0 
is a part of the Schwarzschild space-time. 

Depending on the distribution m(V) of matter in the 
null cloud, either black holes or naked singularities may 
develop as the classical final state of the collapse. For 
example, in the self similar model in which the mass 
is a linear function of the advanced null coordinate, 
m(V) = XV, one finds that both outcomes described by 
the Penrose diagrams in figure 1 are possible, depending 
on whether A > 1/16 (black hole) or A < 1/16 (naked 
singularity) (llj. When the collapse evolves toward a 
naked singularity, spatial hypersurfaces in the future of 
the initial singularity cross the Cauchy horizon and col¬ 
lide with the central singularity but, because no sensible 
boundary conditions can be specified on a singularity, 
the evolution in the future of the initial singularity is ar¬ 
bitrary. The Cosmic Censor should come into play 
before the Cauchy horizon has a chance to form. It is 




FIG. 1: Black hole (left) or naked singularity (right) forma¬ 
tion in null dust collapse 


of interest, therefore, to understand how the system be¬ 
haves close to, but in the past of, the putative Cauchy 
horizon, where spatial hypersurfaces are well defined and 
the quantum evolution of the system may be studied. 

Null shells classically collapse to form covered singu¬ 
larities, therefore, in order to examine such issues as the 
Cosmic Censor, mass distributions other than those rep¬ 
resenting shells should be considered. These require a 
quantization of a genuine field theory. The present paper 
is intended as a first step in this direction. The solution 
metric is written in Eddington-Finkelstein coordinates 
and the Kuchar transformation from ADM variables to 
Kuchar variables is established. The mass function is ex¬ 
plicitly related to m(V ) and m(U) of the metrics in (p~j|). 
respectively s 

The quantization program employed in this paper in¬ 
volves a gauge fixing and it is known that quantum the¬ 
ories resulting from different gauge fixings are not neces¬ 
sarily equivalent. It is, nevertheless, the best approach 
to the quantization of realistic and pressing problems 
such as gravitational collapse at this time. Our choice 
of configuration space coordinates presents several ad¬ 
vantages: their physical meaning is transparent, we are 
able to give explicit transformations between the origi¬ 
nal ADM phase space and the new, and the constraints 
expressed in terms of the new phase space variables are 
linear. The last enables us to obtain solutions for null 
dust clouds of arbitrary mass distributions. The second 
implies that operators representing observables, if known 
in one system are easily constructed in the other, and the 
first makes our solutions easily interpretable in physical 
terms. 


In section |n] we summarize the canonical formulation 
of the action in ADM variables. The null dust action 
appropriate to the models being considered is also ana¬ 
lyzed in this section. In section III we explicitly perform 
a transformation of the phase-space to Kuchar variables 
]l2| . In section IV we apply Dirac’s quantization pro¬ 
gram to these models. Solutions to the constraints for 
collapsing and expanding clouds are presented and the 
matter distribution representing a single shell is exam¬ 
ined as a special case. However, the phase-space admits 
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non-trivial boundaries and this generally leads to compli¬ 
cations in the quantization procedure. Our solutions are 
valid subject to the condition that a suitable measure can 
be found so that the replacement Px —> — iS/5X (r) leads 
to self-adjoint operators. We conclude in section ^ with 
a discussion of the strengths and weaknesses of our ap¬ 
proach, the issues that remain to be resolved (such as the 
question of observables and of the measure on the Hilbert 
space) as well as some suggestions for future directions. 


II. CANONICAL FORMULATION IN ADM 
VARIABLES 


Consider the line element, da on a spherically sym¬ 
metric three dimensional Riemann surface, £. It is com¬ 
pletely characterized by two functions, L{r) and R(r) of 
the radial label coordinate 


dsf 3) = L 2 {r)dr 2 + R 2 (r)dn 2 , (2.1) 

where O is the solid angle. The angular coordinates play 
no role and will be integrated over. We take both L(r ) 
and R{r) to be positive definite except, possibly, at the 
center. R(r) represents the physical radius of a shell la¬ 
beled by r on the surface. It behaves as a scalar under 
transformations of r, whereas L(r) behaves as a scalar 
density. The corresponding four dimensional line element 
may be written in terms of two additional functions, the 
lapse, V(t,r), and the shift, N r (t,r), as 

ds 2 = N 2 dt 2 — L 2 (dr + N r dt) 2 — R 2 dfl 2 . (2.2) 


In this spherically symmetric space-time, we will consider 
the Einstein-Dust system described by the action 


S= - 


1 

16-7T 

1 

87T 


d^xyf—g 1Z 


j d 4 x^g e(x)g a0 U a U f3 (2.3) 


where 1Z is the scalar curvature. As is well known, the 
gravitational part of this action can be cast into the form 


S 9 = 


dt 


dr 


P l L + P r R - NH 9 - N r H 9 r 


+ 


(2.4) 


with the momenta conjugate to L and R respectively 
given by 


Pl 


Pr 


—R + N r R' 

-LR- LR+(N r LR)' 


(2.5) 


and where the overdot and the prime refer respectively 
to partial derivatives with respect to the label time, t , 
and coordinate, r. The lapse, shift and phase-space vari¬ 
ables are required to be continuous functions of the la¬ 
bel coordinates. The boundary action, <Sg S , is required 


to cancel unwanted boundary terms in the hypersurface 
action, ensuring that the hypersurface evolution is not 
frozen on the frontiers. It is determined after fall-off con¬ 
ditions appropriate to the models under consideration are 
specified. The super-Hamiltonian and super-momentum 
constraints are given by 

P L P R LPl 
R 2 R 2 

L R' 2 /RR'\ r 

~2~H + ) 

n 9 r = + R'Pr - LP' l . (2.6) 

We will assume that the matter distribution is such that 
at infinity Kuchaf’s fall-off conditions JiJ| are suitable 
and we will adopt them here. These conditions would be 
applicable, for example, in models in which the collapsing 
metric asymptotically approaches or is smoothly matched 
to an exterior Schwarzschild background at some bound¬ 
ary. They read 


L(t, r) 

= 1 + M+ity- 1 + 

R(t,r) 

= r + 0°°(r~ e ) 

PL(t,r) 

= 0°°(r~ e ) 

PR(t,r) 

= C* 00 ^- 1 " 6 ) 

N(t,r) 

= N+(t) + 0°°(r ~ e ) 

N r {t,r) 

= 0°°(r~ e ) 



and imply that the asymptotic regions are flat with the 
spatial hypersurfaces asymptotic to surfaces of constant 
Minkowski time. Again, as r —> 0 we require that 0 


L(t, r) = L 0 (t) + 0(r 2 ) 

R(t,r) = Ri{t)r + 0(r 3 ) 

PL(t,r) = P L2 (t)r 2 + 0(r 4 ) 

PR{t,r) = PrA^t + 0(r 3 ) 

N(t,r) = N 0 (t) + 0(r 2 ) 

N r (t, r) = N[(t)r + 0{r 3 ) (2.8) 


With these conditions, it is easy to see that the appropri¬ 
ate choice of surface action involves only the contribution, 


S 9 as = ~ [ dtN + (t)M+(t) (2.9) 


at the boundary at infinity. 

Let us now consider the null dust action in (2.3). We 
note first that the energy density, e(x), plays the role of 
a Lagrange multiplier enforcing null dust, i.e., U 2 = 0, 
and variation w.r.t. gap yields the standard dust stress 
tensor, T a p = e(x)U a U^. The canonical form of the null 
dust action in various forms has been studied by Kuchar 
and Bica k |fl5[ . In particular, for the action in the form 
given in ( p.3| ) one may expand U a as a Pfaff form of six 
scalar fields, the three co-moving coordinates of the null 
dust particles, Z k , and three scalars (velocities), Wk, 


U a = w k Z k 


( 2 . 10 ) 


This representation is redundant because, by Pfaff’s the¬ 
orem, only four scalars are required to describe an ar¬ 
bitrary covector in a four dimensional space. Suppose 
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we require one of the scalars, say W 3 to be unity and 
drop the index from the associated co-moving coordinate, 
Z 3 := Z , then 

U a = Z, a + w k Z k a , k G {1, 2}. (2.11) 

Consider the independent variations, 


Setting SC/Se = 0 gives the final form of the dust Hamil¬ 
tonian and momentum constraints 


n d 

njd 

ri r 


L P z (Z'+ w k Z k ') 
L 

P z (Z'+ w k Z k '), 


(2.18) 


° = -g = -V=gg al3 u a u 0 

0 = = V a (e(x)U a ) 

0 = 

° = = (V^ge{x)w k U a ), a 

= V a (e(x)w k U a ) 


where the positive (respectively negative) signs in the 
dust Hamiltonian density represent incoming (respec¬ 
tively outgoing) dust. In the spherically symmetric col¬ 
lapse we are considering, we take w k = 0 = P k . Thus we 
have arrived at the canonical form of our theory, which 
we will write as 

S = + j dt j dr[ZPz + LPl+ 

+ RP r - NH - N r Hr\ 


2 6 S 

T a/ 3 = — 7 - = e(x)U a Up (2-12) 

V~9 °9a0 

The conservation of the stress energy tensor in the last 
equation implies that 

e{x)U a W a U 0 + U 0 X7 a (e{x)U a ) = 0, (2.13) 


which says that the particles follow geodesic curves. Us¬ 
ing the second equation above we find U a \7 a U@ = 0, im¬ 
plying affine parameterization. The third equation says 
that CjjZ k = Z k a U a = 0, i.e., all of the Z k are constant 
along flow lines and none of them are time-like. And 
finally, multiplying the third equation by w k we find 

Z, a U a = 0, (2.14) 


saying that Z a may be space-like or null. If the twist, 
U[ a .f 3 ], also vanishes, then Z >a is null, which would imply 
that w k Z k a = 0, ov w k = 0 V fcG {1,2}, because the Z k a 
are taken to form a (linearly independent) cobasis. 

Substituting the decomposition (2.11) into the dust ac¬ 
tion in (2.3), using (|T^) and integrating over the angular 
coordinates the action may be put in the form 


S d 



P z Z + P k Z k - NH d - N r U d r , 

(2.15) 


where the momenta conjugate to {Z, Z k } are, respec¬ 
tively, 


Pz = -^-e(r)[(Z + w k Z k ) — N r (Z' + w k Z k ')} 

Pk = w k Pz, (2-16) 


and the constraints, Ti d and H d are 


= \ P 2 z , eR\Z'+w k Z k 'f 
2 LR 2 e 2 L 

U d = P z (Z'+ w k Z k '). (2.17) 


H = 


n 9 r = 
Sas = 


PlPr LP 2 l 
R 2 K 2 

L R ' 2 (RR'\ 

~2 ~ ~2L + {—J 


Z'Pz + R'Pr - LP' l 
[ N + (t)M + (t ), 


Id So 


- V 


PzZ' 

L 


(2.19) 


where 77 = sgn(Z r ). In the following section we will show 
that the co-moving coordinate Z may be identified with 
the null coordinates according to Z = — U for an expand¬ 
ing solution and Z = —V for a collapsing one. Pz < 0 
and the dust Hamiltonian density is chosen to be always 
non-negative. When Pz is non-vanishing the phase-space 
is made up of two disconnected sectors, labeled by 77 . An 
initial data set with Pz = 0 cannot evolve into a set with 
Pz 7 ^ 0 and we will assume from now on that Pz 7 ^ 0 . 


III. CANONICAL TRANSFORMATION 

The description of contracting and expanding clouds is 
seen to be related by time reversal. The two descriptions 
may be formally unified in the following way. Introduce 
a null coordinate W n . which can be the “advanced” time 
or the “retarded” time, satisfying only the requirement 
that W v increases toward the future. If W' v < 0 (primes 
denote differentiation w.r.t. the ADM label coordinate 
r) then it represents the retarded coordinate, U , and if, 
on the contrary, W' > 0, it represents the advanced co¬ 
ordinate, V. Let us write both solutions in terms of a 
parameter 77 that represents the behavior of the matter 
(whether it is expanding or collapsing) 

1 - — j dW r 2 + 2-ridW ri d R - R 2 C m 2 . (3.1) 

The metric is appropriate for either expansion or 

contraction of the dust cloud depending on whether 77 = 
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—sgn(W,/) is +1 or — 1 . rj is the same as appears in 
Q , as we argue below. The null coordinate, W v , 
whose spatial direction is opposite to W r/ is obtained by 
integrating 


dW v = a(W v ,R) 


dW v + 2 r] 


dR 

~T 


(3.2) 


where A = ( RR')(LPl )' — (RR')'(LPl). One can then 
show that the transformation, 

(R, L, Pr, Pl) —> (R,m, Pr, Pm), (3-11) 

is canonical, and generated by 


where a(W v ,R) is an integrating factor. This coordinate 
must also be always i ncre asing toward the future. 

The hypersurfaces dU from which ( |2.2| ) is constructed 
must be em bedded in the space-time described by the 
metric (3.1}. Substituting the foliation W v (t,r) and 


R(t,r) in (|3.l|) gives the density L(t,r ) and the lapse 


and shift functions, N(t,r ) and N r (t,r), as 


TWl + 2 V W v R 
-TWrJ 2 - 2riW r /R' 
-V{W V R' + W V 'R) - FW V W V ' 


N 2 - L 2 N r2 
L 2 

N r L 2 , (3.3) 


where we have set T = 1 — 2m/R. These relations can 
be used to determine, 


PW V W V ' + v(W v R' + W V 'R) 
PW V ' 2 + 2 r/Wr/R! 
W V R' - Wr/R 
L ’ 


(3.4) 


where we have chosen the sign of the square root so that 
N is positive. Inserting these expressions for N and N r 
in the expression ( p.5| ) for Pl, we find 

^ = -VR' - TWr !, (3-5) 


whic h, when substituted into the expression for L 2 in 

(0, gives 

/2 


T = 


R! A 


Pl 


L 2 R 2 ’ 


(3.6) 


or, equivalently, the mass function in terms of the canon¬ 
ical data 


G 



LPl 


X -RR' In 


RR' + LP l 
HR' - LP l \ 


(3.12) 


By computing the difference between the old and the 
new Liou ville forms and using the fall off conditions in 
© and ( | 2 . 8 | ), one can show that the transformation has 
introduced no fresh boundary terms. 

There are (infinite) boundary terms at the horizon, 
when T — 0. It can be shown, however, that the con¬ 
tribution from the interior and the exterior cancel each 
other. There will also be contributions at the boundary 
between the interior of the star and its exterior or more 
generally at any frontier between two regions described 
by mass functions with different derivatives. Again, if 
the mass function is continuous across the boundary and 
regions are consistently matched by equating both the 
first and second fundamental forms, then the contribu¬ 
tion from one side will cancel the contribution from the 
other. 

Before rearranging the action, we will consider the co¬ 
ordinates {Z,Z k } for the collapsing Vaidya null congru¬ 
ence. Returning to the metric in 0 ) with 77 = — 1 
we find that the incoming null congruence is given by 
V = const., 6 = const, and <f> = const. The coordinates 
Z = —V, Z 1 = 9 and Z 2 = <j> are co-moving. Let us form 
the basis 


z,n = (- 1 , 0 , 0 , 0 ) 

Z^ = ( 0 , 0 , 1 , 0 ) 

Zl = (0,0,0,1) (3.13) 



ir_ 

H 2 



(3.7) 


By directly taking Poisson brackets, the momentum con¬ 
jugate to m can now be shown to be simply 


Prr 


LPl 
RT ' 


(3.8) 


Kuchaf |jl 2 fe proposed that (R,m, Pr, Pm) should form 
a canonical chart whose coordinates are spatial scalars, 
whose momenta are scalar densities and which is such 
that H r (r) generates Diff R. This means that 


LLr = P'Pr ^ LP/ = R'Pr + m P m ~ 0 . (3.9) 

Substituting the expressions derived for to and P m into 
the above constraint one arrives at 


LP l LP l A 


(3.10) 


It is easily shown that £ = — R is an affine parameter 
and the covariant components of the velocity dx^/d^ 
are = (— 1 , 0 , 0 , 0 ), whose decomposition in the co¬ 
basis Z k a yields W\ = 0 = Wi- Similarly treating the 
outgoing null congruence shows that the affine param¬ 
eter is £ = +R and that Z is to be identified with 
—U. Both cases may be treated simultaneously by letting 
Z = —W v , Pz = —Pwr, an d £ = tjR- This identification 
shows that the 77 used in the section is identical to that 
used in the previous section. 

N ote also that taking the spatial derivative of to in 
( 0 ), yields 


m ' = ~T n 9 -^ n r- 


(3.14) 


This may be used to write the action in (2.3) as 


Pr = Pr 


2 R 2 RT RL 2 T ’ 
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S = J dt J dr 


Pwn W v + PrR + P m rn - NTL - N r H r 

I- 


SdY ., 


(3.15) 


with 


n = - 


m'T^R! + TP m P R 


Pw v W v ' 


i.e., m'(r) = F^(r) and M + is clearly the mass at infinity. 
The one form 


7 H r — R' P r + rn'Pm + W v 'P Wv 
SflE = - / N+(t)M+(t). 

JOE™ 


(3.16) 


fl= drP m (r)Sm(r) — M + SW V+ 

Jo 


can now be written as 


(3.19) 


The surface term contains the mass at spatial infinity and 
may be re-expressed in a more convenient form. Use the 
fall-off conditions (2.7) at infinity and the expression for 
N in (3.4) to write N + = W v+ and 


= - / W v+ (t)M + (t). (3.17) 


IdS o 


Then define T r) (r) by 


n = + J drP m (r)J SM + 

POO POO 

— / drP m (r ) / dr'Sr v (r') — 6u , (3.20) 

Jo J r 

where 6ui = S(M + W V , ) is an exact form. The second 
he expressioi 
be cast into 

identity 


term in the expression lor 0 continues to be inconvenient, 
m(r) = M + — / drl'r)(r ! ) (3.18) but may be cast into a more appropriate form using the 

J r 


pr p oo \ f poo / r r \ 

J ' dr'P m (r') x J dr'5T v {r' )j = P m (r) J dr'5T v (r') - dr’P m {r')UT v (r). (3.21) 


r 


Integrating from r = 0 to r = oo, the left hand side of where we have defined 
the above equation vanishes identically and one finds 


p OO POO 

J dr P m (r) J dr'ST v (r') = 

drSTr,(r) f dr'P m (r'),( 3.22) 
Jo 


P+ = W v+ + j drP m (r ) 


jo Jo 

so that fl can be cast into the form 


Pr n = ~ I dr'P m (r'). 
Jo 


(3.24) 


POO 

= p+5M + + / drPr ST^ — Sui, (3.23) Eliminating a total time derivative turns the action in 

Jo v ( |3.15 ) into 


S = dt( p+M+ + / dr[P w W ri + P r R + P Tv r, -NH- N r H r ] , 


(3.25) 


where 


n = - 


r^- 1 


R'-PP'P f 


w. 
























7 


Hr = r'Pr - r ri p^ + w r ;p Wv . 

Furthermore, it follows from (|T^) that 

R' 


P^=W rl ' + T 1 


T 


(3.26) 


(3.27) 


where <j(R, W rl ) is the integrating factor introduced in 

the previous section. But W v ' ^ 0 because it is a null 
coordinate and is required to increase in time, therefore, 


Ph-V^^O 


(3.30) 


The constraints, Ti sa 0 ~ 7i r , can be further simplified 
by using H r ss 0 to eliminate W v 'Pw n from the Hamilto¬ 
nian constraint. This gives 



(3.28) 


Consider the first of the two factors above. Using ( 3.27 ) 
to substitute for Pp , we find 


TP^ + r,m = T ( W r ; + 2^) = 


(3.29) 


is equivalent to the Hamiltonian constraint, Ti ~ 0. In¬ 
serting this into either of the two constraints then gives 

Pw n ~ v ppr ~ r,,. 

The configuration space is made of the set of variables 
{ W rj . R, r v , M + }, whose physical significance is transpar¬ 
ent. This is an advantage of Kuchar variables. W v is a 
null coordinate, R is the area radius of a point labeled 
(r,t), T,, is the energy density of the collapsing cloud 
and M + is the mass measured at spatial infinity. M + 
is a constant of the motion and may be viewed as part 
of the initial data. Our gravity-matter system may be 
re-written as 


J 


/ 


S = dt[ p+M+ + 


dr 


Pw„W n + PrR + Pr^ri 


The canonical action (3.31) can be a starting point for 
Dirac quantization. The physical meaning of all the con¬ 
figuration space coordinates is clear: W v and R locate 
the hypersurface and T^ (along with M + ) determines the 
matter distribution. 


IV. QUANTIZATION 

The configuration space consists of two disconnected 
components, for expanding null matter it is spanned by 
the set {W + , R,T + M + } and for collapsing matter by the 
set {W-, R, T_, M + }. In each case the phase-space has 
non-trivial boundaries, 

M+ >0, P z < 0, (Z = -W±) 

R > 0. (4.1) 

The last is due to the fact that R(t, r) = 0 describes the 
central singularity, i.e., the final singularity for the col¬ 
lapsing solution and the initial singularity for expanding 
one. In Dirac’s approach, when the phase-space admits 
trivial boundaries, the canonical momenta Px , are raised 
to operator status, 


Pv - -4 (4.2) 

and the constraints are considered as operator restric¬ 
tions on the state functional. When the boundaries are 


N w (P Wri - T v ) - N r (R'P r - T„P' + W v 'P w ) 


(3.31) 


r 


non-trivial, as is the case in (4.1), this naiv e exchange of 
momenta for functional derivatives in (4.2) may lead to 
operators that are not self-adjoint [^3|, but we will as¬ 
sume here that a suitable measure on the Hilbert space 
can be found, with respect to which the above replace¬ 
ment leads indeed to self-adjoint operators. Subject to 
this caveat, the state-functional obeys 


AvJ/ Ailr 

R'^Z + w'— - r„ 


SR 


SWr, 


. S 4' 

"Wr, 


(5T 

sr: 


= 


= 0 , 


(4.3) 


The last constraint is solved by any functional that is a 
spatial scalar. Consider a solution of this constraint that 
is of the form 


'F = C rj {M + ) exp 


/»oo 

i / drT v (r) • ICfjiiWrj, R,M+) 

Jo 


(4.4) 


where C v is a constant depending only on 77 and M + , 
and K, is an arbitrary complex valued function of its ar¬ 
guments (and not their derivatives) that is to be evalu¬ 
ated so that T satisfies the other constraints. The wave- 


functional T in (T4) is evidently a spatial scalar because 


T ? 7 (r) is a spatial density and 1C is a spatial scalar. It 
is therefore a solution of the constraint providing that 1C 
has no explicit dependence on r. 

The solution, which agrees with all the constraints con¬ 
sidered as operator restrictions on the state functional, 
















is given by 


4/ = C v exp 


drT v (r){W v (r) +r]R*(r)) 


(4.5) 


where R^,(R, M + ,T V ) is a “tortoise”-like coordinate de¬ 
fined by 


R 

R* = R + 2 m In I- 1 1. 

2 m 


(4.6) 


It is not, of course, the tortoise coordinate R* except in a 
Schwarzschild region when m(M + , T,,) = M + is constant. 

The parameter 77 represents the direction of the flow, 
being +1 for outgoing null matter and —1 for incoming 
matter. The combination r]R(r) represents the affin e pa - 
rameter, £(r). Re-expressing the wave-functional in (4.4) 
to make the dependence on the affine parameter explicit, 
we find 


T = C v e^ exp 


drT v (r) ■ W v (r) 


(4.7) 


where £* = r]R*. 

For 77 = — 1 , £ £ (— 00 , 0) and W v = V, the functionals 
are of the form 


^r _ 1 = A(M + )e~^o° drr -( r ) R *( r ) e i f™ drT—(r)V(r) 


(4.8) 


and describe collapsing null matter. Likewise, for 77 = 
+1, £ € (0, 00 ) and W v = U, the functionals 

4»+i = B{M + y5™ drT +( r)R * {r) e i 5~ drT+{r)u(r \ (4.9) 


describe expanding null matter. 

A given classical collapse problem is specified by a 
choice of mass function, m(W T? ), which determines an 
initial energy distribution, thus a collapse “model”. As 
an example, we shall consider the function, 

m(W v ) = M + e{W v - w), (4.10) 

where 8 is the Heaviside unit step-function and w is con¬ 
stant. The matter energy vanishes when < w and is 
M + when > w. The mass function evidently makes 
sense only as a thin shell that is collapsing toward the 
center and we must have W v = V (77 = —1). Calling the 
corresponding mass at spatial infinity MU, we find the 
energy density by differentiating w.r.t. the ADM label 
coordinate, r, 

r mn (r) = M?V'(r)6(V-v) 

■ V'(r) 

= (4.11) 

where t(t) is the solution of V(r, t ) = v. Likewise, a thin 
shell that expands out of the center is represented by the 
mass function 


m(W v ) = Ml ut 8(w - W v ), (4.12) 



for W v = U (77 = +1). The energy density is 


^out 


(r) = -M“ ut l/'(r)(5(« — U) 


TJ'(r) 

= M™» t 777 Ll ( J(r-t), 


U'( t) 


(4.13) 


where r (t) is the solution of U(r,t) = u. The energy 
density is always positive and in either case we find 


T = M+ n(out) <5(r - r). 


(4.14) 


The shell trajectories, r = r (t) are different in the two 
cases, being v = const, in the first and u = const, in the 
second. Us ing t hese expressions, let us relate the con¬ 
straints in ( [2.19 ) to the constraints that have been used 
by others to describe thin shells. The dust Hamiltonian 
and momentum density turn out to be (Pz = — T) 

TL d = 77 — j — S(r — r) = r/yS(r — r) 

JL/ Lj 

= —M + Z'S(r - r) = -p S(r - r), (4.15) 

where we have defined p = M + Z'(x) and therefore 77 = 
sgn (Z') = sgn(p). These expressions were used as a start¬ 
ing point in ;[T~l| . 00 - The gravitational contributions 
to the constraints are, of course, the same. 

The corresponding classical solutions are represented 
in the Penrose diagram of figure 2, where A A' represents 
the event horiz on in ( a) a nd the Ca uchy horizon in (b). 
Inserting ( 4.14 ) into (4.8) and (|4.9| ) one finds that the 
quantum mechanics of a single collapsing shell is defined 
on the reduced configuration space (v,R,M+) and de¬ 
scribed by the wave-function 


$_! = A(M; n ,u) e - lM i" fl * (t) , 


(4.16) 


whereas, for an outgoing shell the reduced configuration 
space is {u,R, M? ut ) and the wave-function 


tf +1 = 5(tf” ut ,«)e iM ” ,fl *W, 


(4.17) 


where R* as given by (4.6) is, in this case, the usual 
tortoise coordinate. 
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V. DISCUSSION 


In this paper we have examined the collapse and ex¬ 
pansion of a null dust cloud of arbitrary mass distribu¬ 
tion and shown that there exists a canonical transfor¬ 
mation that brings the corresponding Vaidya system to 
the Kuchar form, in which the dynamics is expressed in 
terms of embedding variables whose physical meaning is 
transparent. Written in terms of these variables, the con¬ 
straints take on a simple form. 

The physical content of the solution wave-functionals 
discussed in the previous section is given by the expec¬ 
tation values of observables. These are, for example, the 
geometric invariants and matter invariants (such as cur¬ 
vature scalars and trace of the stress energy tensor) and 
in general any C°° function of the phase-space variables, 
old or new, that weakly commutes with the constraints. 
They must be written as operators on the Hilbert space, 
but there are two associated difficulties. Firstly, because 
they are generally non-linear, there will be operator or¬ 
dering ambiguities. Secondly, one must ensure that they 
are self-adjoint with respect to the chosen measure. How¬ 
ever, as our transformation between the spaces is explicit, 
knowing these functions in one system is equivalent to 
knowing them in the other. It should be noted that a 
complete set of Dirac observables has been constructed 
in the case of a single shell [^6| where a single degree 
of freedom is present. For the general case, this issue is 
quite complicated and will be discussed in the future. 

Subject to the condition that suitable boundary con¬ 
ditions or a suitable inner product can be found so that 
(4.2) leads to self-adjoint operators, we have found so¬ 
lution wave-functionals for arbitrary mass distributions 
in each sector (collapse and expansion) independently 
(as mentioned in the introduction, solutions with non¬ 
trivial mass distributions are essential for the description 
of issues as naked singularities and the Cosmic Censor). 
These sectors are disjoint, separated by the central sin¬ 
gularity (at R = 0). However, if the solutions are viewed 
as describing the evolution, in the affine parameter £, of 
a collapsing matter cloud beginning on 3", they seem 
to suggest that it should be possible to define the wave- 
functional over the entire interval £ G (— 00 , 00 ) by ex¬ 
tending the range of R to include the center. Yet , fo r 
any matter distribution, the solution space-time in ( |3. l| ) 
admits a strong curvature singularity at the origin, so the 


classical dynamics cannot be extended even to it. Any 
attempt to continue the quantum dynamics through the 
origin must therefore ensure that at least the expecta¬ 
tion values of the observables in the consequent quantum 
theory are well behaved there. Thus the central singu¬ 
larity would be made harmless by the quantum theory. 
It can be so if, for example, the wave-functional were 
to vanish there. In this way the matter would collapse 
and re-expand through the (benign) center in one con¬ 
tinuous history, the solutions being given by ([h^) for 
£ € (—oo,0] and (4.S) for £ £ [0, 00 ). This has been pro¬ 
posed for a single shell by Hajicek and Kiefer El EMI. 
who merged the two solutions into one bouncing solution. 
Their bounce was obtained by working in double-null co¬ 
ordinates and employing group quantization techniques 
(see |nj and ) to the problem. Group quantization is 
beautifully adapted to the quantization of systems with 
non-trivial boundary restrictions such as those in (4.1) 
on the phase space, but it’s application to problems with 
more than a few number of degrees of freedom and in 
particular to the collapse of general matter distributions, 
being dependent on the construction of a complete set of 
observables, remains a topic for future investigation. 

The Eddington-Finkelstein coordinates we have em¬ 
ployed in this paper present several advantages over the 
double null system in regard to the problem of collapse or 
re-expansion of arbitrary matter distributions. The new 
variables have a clear physical and geometrical mean¬ 
ing. This is useful when comparing the quantum behav¬ 
ior with the classical. Our transformations are explicit, 
which means that operators and, in particular, observ¬ 
ables that are known in one coordinate system can be 
expressed in the other. The matter-gravity constraints 
in the new phase space are linear for all matter distribu¬ 
tions. This simplification, achieved on the classical level, 
has allowed us to obtain exact solution wave-functionals 
describing the respective physical processes. 
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